On the decrease of intermittency in decaying rotating turbulence 
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The scaling of the longitudinal velocity structure functions, Sq(r) — {\5u{r)\'') ~ r""', is ana- 
lyzed up to order g = 8 in a decaying rotating turbulence experiment from a large Particle Image 
Velocimetry (PIV) dataset. The exponent of the second-order structure function, ("2, increases 
throughout the self-similar decay regime, up to the Ekman time scale. The normalized higher-order 
exponents, are close to those of the intermittent non-rotating case at small times, but show a 

marked departure at larger times, on a time scale {Q is the rotation rate), although a strictly 
non- intermittent linear law C<j/C2 = ?/2 is not reached. 



Whether intermittency of isotropic three-dimensional 
(3D) turbulence is decreased or even suppressed in 
the presence of s yste m rotation has recently received 
a marked interest .1321 Here, intermittency refers to the 
anomalous scaling of the structure functions (SF) of 
order q, Sq{r) = (|i5M(r)|'') ^ r*"', where 5u{x,r) = 
[u(x -f r) — u(x)] • r/r is the longitudinal velocity in- 
crement, r an inertial separation normal to the rotation 
vector n and (•) denotes spatial and ensemble average. 
A linear variation of the exponents Cq with the order 
q is the signature of self-similar (non-intermittent) ve- 
locity fluctuations, a situation which is found in the in- 
verse cascade of two-dimensional (2D) turbulence.^ On 
the other hand, anomalous exponents, C0 /C2 7 ^ 9/2, are 
the landmark of 3D isotropic 

turbulence .HiSsI Based on 
the qualitative ground that rotating turbulence experi- 
ences a partial two-dimensionalization, one may naively 
expect a reduction or a suppresion of intermittency by 
comparison with the 3D non-rotating case. More pre- 
cisely, describing rapidly rotating turbulence in the limit 
of zero Rossby numbers as a sum of weakly interacting 
random inertial waves, the vanishing of non-linear effects 
should lea d to a special case of non-intermittent wave 
turbulence. 

Two papers have recently addressed the issue of the 
scaling of the SF in rotating turbulence with a station- 
ary forcing. The hot-wire measurements of Baroud et al.^' 
in a turbulent flow generated by radial jets in a rotating 
tank showed a transition from an intermittent to a non- 
intermittent behavior, characterized by a E{k) ~ fc~^ en- 
ergy spectrum (i.e. C2 = 1) and linear higher-order expo- 
nents (^q — q/2. In a Direct Numerical Simulation (DNS) 
of rotating turbulence with a large-scale isotropic forc- 
ing, Miiller and Thielc3 have observed reduced intermit- 
tency, also characterized with C2 — 1, but higher-order 
exponents (^g intermediate between q/2 and the values 
usually found in classical (intermittent) 3D turbulence. 
Those observations are in qualitative agreement with the 
increase of reported by SimancP' from hot-wire mea- 
surements in the vicinity of a strong vortex, although no 
clear separation between a constant background rotation 
and an otherwise homogeneous turbulence advected by 
the rotation can be defined in this geometry. To date, no 
theoretical description of the scaling of the anisotropic 



higher order SF in rotating turbulence is available. Note 
that, in all the above references, only separations r nor- 
mal to the rotation vector ft are considered, ignoring the 
complexity originating from the anisotropic character of 
rotating turbulence. 1^ 

In this Letter we report new measurements of the 
high order SF, carried out by Particle Image Velocimetry 
(PIV), in a freely decaying rotating turbulence experi- 
ment, aiming to compare with the results obtained in 
forced turbulence. The experimental setup is the same 
as in Morize et al.,'^ and is only briefly described here. 
It consists in a water filled glass tank of square section, 
of side 35 cm and height /i = 44 cm, rotating at constant 
angular velocity. After the fluid is set in solid body ro- 
tation, turbulence is generated by towing a co-rotating 
square grid, of mesh size M = 3.9 cm, at a constant ve- 
locity Vg — 0.65 m s~^ from the bottom to the top of the 
tank, and is maintained fixed near the top during the 
decay of turbulence. The horizontal components of the 
velocity fields in a centered horizontal area of 17 cm x 14 
cm at mid height of the tank are obtained using a coro- 
tating PIV system operating at 1 Hz. The velocity fields 
are defined on a 160 x 128 grid, with a spatial resolution 
of 1 mm and a signal-to-noise ratio of about 2 x 10~^. 
Although this fails to resolve the dissipative scales (the 
Kolmogorov scale is approximately 0.2 mm in the first 
period of the decay), this resolution allows us to resolve 
the inertial range, typically for r > 10 mm. 

Two angular velocities have been used in the present 
experiments, = 1.13 and 2.26 rad.s"^. The correspond- 
ing non-dimensional parameters are summarized in Ta- 
ble |l] The grid Reynolds number is Rcg = VgM/v — 
2.5 X 10'' {v is the kinematic viscosity), and the grid 
Rossby numbers Rog — Vg/{2Q.M) are 7.4 and 3.7, so 
that the initial state can be considered as a fully devel- 
oped 3D turbulence weakly affected by the system rota- 
tion. A previous investigatioiP^I showed that, for those 
rotation rates, the energy decay was approximately self- 
similar between tp ~ 40M/T^ and — 0.10^^;, where 
Ie ~ h/{v^Y/'^ is the Ekman time, followed by an expo- 
nential decay at larger times. The present investigation 
is restricted to the self-similar range [Iq , tc] ■ The instan- 
taneous Reynolds, macro- and micro- Rossby numbers. 
Re M , Rom and Ro^^ respectively, are also given for these 
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FIG. 1: Time evolution of the second order exponent for ^ = 
1.13 rad s~^. — , whole velocity field; o, turbulent field (ensemble 
average substracted). 



two limiting values to and tc in Table [Tj 

To ensure proper convergence of the statistics, each 
decay is repeated approximately 600 times, represent- 
ing 10 hours of run for each rotation rate. It is worth 
pointing that computing SF from PIV data requires spe- 
cial care, especially when higher order are considered, for 
which even a small number of spurious vectors may have 
a large effect. Since those bad vectors may be preferen- 
tially found in regions of large velocity or large gradient, 
finding correct criteria for removing them without intro- 
ducing biases is a delicate issue. In particular, some of 
the fields were found to suffer from an inhomogeneous 
lighting, because the imaged area was partially shad- 
owed when the corner of the tank passed through the 
laser sheet. Using a criteria based on the Q-factor (ratio 
of primary and secondary correlation peak) , 20% of the 
fields were affected by this problem and have been re- 
moved. A median filter is then applied to the remaining 
fields, and it was checked that the SFs computed from 
the raw and median-filtered data agreed within the error 
bars ASq defined below. 

We first focus on the time evolution of the exponent 
C2 of the second order SF, S2{r) = {\Su{r)\^), plotted in 
Fig. [1] for — 1.13 rad s^^. This exponent is related 



TABLE I: Non-dimensional parameters for the two rotation 
rates, f tp ,tc] is the range of approximately self-silimar energy 
decay.E21 Rsg and Rog are the grid Reynolds and Rossby num- 
bers. ReM{t) = u'{t)M/u, RoM(t) = u'{t)/{2Q.M) and Ro^(t) = 
uj'{t)/2Q are the instantaneous Reynolds, macro- and micro- 
Rossby numbers respectively, based on the horizontal velocity rms 
u'(t) and vertical vorticity rms Lo'it). 



n (rad s-i) 



mo/27r . . . ntc/2TT 
ReM{t = U) ■ ■ ■ tc) 
RoM{t = to . . .tc) 
Ro^{t = to...tc) 



1.13 

2.5 X 10"' 
7.4 

1.2 ...7.4 
1300 ... 360 
0.38 ...0.10 
2.1 . . .0.23 



2.26 

2.5 X lO"* 
3.7 
0.6 . 
1400 
0.21 
1.1 . 



10.5 
. .380 
. 0.056 
0.17 



to the distribution of energy among scales: larger val- 
ues of C2 indicate a favored energy distribution toward 
larger scales. Significant oscillations of (,2 are present, 
with a period equal to the tank rotation period, indi- 
cating the presence of inertial modes. Those inertial 
modes have been previously detected from oscillations 
in the decay of the kinetic energy by Morize et al.,'^^! 
and their temporal spectrum has been analyzed in de- 
tails by Bewley et al!^ Since we are interested here in 
the turbulent fluctuations that superimpose to those slow 
modes, we have computed the turbulent velocity fields 
u"(x, t) = u"(x, t) — (u"(x, t))a, where a denotes the re- 
alization and (•)q. is the ensemble average over the whole 
data set at a given time t after the grid translation. The 
time evolution of the corrected exponent C21 measured 
from the scaling of the turbulent component of the SF, 
S2{r) = {\Su{r)f) (also plotted in Fig. [ij, is found to 
follow approximately the lower bound of the oscillations 
of the raw exponent (^2 ■ One may conclude that the iner- 
tial mode, by superimposing a large scale modulation to 
the turbulence, leads to an increased raw exponent (,2, of 
order of 10%. In the following we will discard this slow 
inertial component of the fiow and we will focus on the 
scaling of the turbulent flow component. 

The corrected exponent, hereafter simply noted C2, is 
found to gradually increase during the decay, starting 
from values close to 2/3 at t ~ to, as expected for an 
initial state weakly affected by rotation, and increasing 
up to 1.4 ± 0.05 at t ~ tc, reflecting the growing im- 
portance of the large scales compared to the small ones. 
This behavior compares well with the gradual steepening 
of the energy spectrum reported by Morize et al.,'^with 
a spectral exponent p increasing from 1.7 to 2.3 ±0.1 
during the decay [dimensional analysis gives C2 = P ^ 1 1 
with E(k) ~ k~P the one-dimensional spectrum com- 
puted from the horizontal velocity and k the horizontal 
wavenumber]. Beyond tc, the energy decreases exponen- 
tially as the result of the dissipation by the inertial waves, 
and no s caling range could be defined from the power 
spectrum .tlillll In the following we restrict to times t < tc, 
where a correct scaling over an appreciable range of scales 
is observed from both S2{r) and E{k). 

We now turn to the higher order SFs. Figure |2ja), 
where SFs up to order q = 8 are plotted at a given 
time t, show power laws for intermediate scales, here for 
12 mm < r < 80 mm. It is worth pointing that the de- 
termination of the highest measurable order and its un- 
certainty for a given sample size is a delicate issue. The 
highest order for converged SF is determined by visual 
inspection of the truncated integral 



p{5u) \5u{r)\'^ dSu, 
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which increases up to Sq{r) as the cutoff Su* 
increased.'^ Here p is the probability density function 
(pdf ) of the velocity increment Su. For large separations 
and/or moderate orders, Cq increases smoothly towards 
a well defined plateau as 6u* 00, indicating a correct 
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FIG. 2: Structure functions for increasing orders, for Q = 
1.13 rad s~^ and Q,t/2-K = 2.7, normalized by the velocity rms 
u'(t) (a) plotted as a function of the separation r; (b) plotted as a 
function of 52 (ESS method). The curves for g = 4, 6, 8 have been 
vertically shifted by factors 10'^, 10* and 10® for visibility. The 
dashed lines show the range where the exponents are fitted. 



convergence of the SF. On the other hand, smaller separa- 
tions, r < 10 mm, show strong jumps when large velocity 
increments enter into the integral ([T|) . Those jumps may 
be due to either spurious vectors or insufficent statistics, 
and are the signature of an unconverged SF. According 
to this criterion, the range of separations r ensuring a 
correct convergence of 5*^ for g > 8 is found too small 
for a reliable measurement of the scaling exponents, and 
measurements are restricted to order q = For orders 
g < 8, scales r > 10 mm were always correctly converged, 
allowing to safely define scaling exponents in the inertial 
range. Finally, the uncertainty AS'q(r) is estimated by 
plotting Sq(r) at a given order and a given separation 
as a function of the sample size. Defining ASq{r) as the 
standard deviation of Sq{r) computed over the last third 
of the whole sample yields a relative error l^Sq{r) / Sq{r) 
of 4% for q = 4 and 10% for q = 8, which is smaller than 
the symbol size in Fig. [2j 

Figure [31 shows both the raw exponents C,q and the 
normalizea exponents Cg/C2 at different times during the 
decay. Those raw (resp. normalized) exponents are ob- 
tained from a linear least-square fit of log Sq versus log r 
(resp. logiS'2), following the Extended Self Similarity'* 
(ESS) procedure [see Fig. |2]^b)].''* The main contribu- 
tion of the error bars for Qq is due to the uncertainty 
on the determination of the SF discussed above, A^^ ~ 
2(A5'q/S'q)/ln(r2/ri), where ri and r2 are the lower and 
upper cutoffs of the scaling range, yielding A^4 ~ 0.05 
and ACs ~ 0.14. At the beginning of the decay, the ef- 
fect of rotation is small and the exponents are indeed 
found very close to classical values for 3D non rotat- 
ing turbulence .'S^ For comparison, the She-Levequ^ for- 
mula is also plotted, showing good agreement up to order 
g = 8, giving confidence on the reliability of our PIV mea- 
surements. At larger times, the normalized exponents 
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FIG. 3: (a) Raw exponents C,q, and (b) normalized exponents C9/C2 
measured using ESS, at various times during the decay, for £7 = 
1.13 rad s'^ . 



increase and become closer to the linear law Cg/C2 — 9/2, 
confirming the intermittency reduction induced by the 
background rotation. It is worth noting that the instan- 
taneous Reynolds number a,i t ~ tc, Rcm — 360 (see 
Table |l]) , together with the correct scaling of the SF at 
that time, ensure that this intermittency reduction is 
not associated to the trivial scaling Q — q (and hence 
Cg/C2 = 9/2) of a smooth velocity field. 

The exponents at the end of the decay are compara- 
ble or even slightly larger than those reported by Miiller 
and Thiele,^ although their macro- Rossby numbers (0.01 
and 0.05) are slightly lower and their Reynolds num- 
ber (2300 and 4000) significantly larger than the present 
ones (note that the non-dimensional numbers here are 
based on the mesh size M, which underestimates the 
true integral scale). It must also be noted that the 
present exponents differ from the strictly linear law q/2 
reported by Baroud et al} for similar Rossby numbers. 
This slight discrepancy may be due to the different forc- 
ing mechanisms: In the present experiment, the initial 
turbulence produced by the grid translation is approxi- 
mately isotropic, and rotation gradually breaks this ini- 
tial isotropy in the course of the decay. In the experiment 
by Baroud et aZ. turbulence is maintained by radial jets 
originating from a circular array of holes, generating a 
strong radial flow deflected by the Coriolis force. This 
forcing scheme is likely to produce an anistropic, par- 
tially two-dimensional flow, even in the absence of rota- 
tion. When rotation is present, this forcing probably re- 
inforces the two-dimensional character of the turbulence, 
resulting in strictly non-intermittent exponents. 

The reduction of intermittency during the decay is best 
appreciated from the intermittency factors 7^ = q/2 — 
Cg/C2 (Figs. [4] and [5|, which vanish for non-intermittent 
fluctuations. Although the scatter is important on these 
quantities (of the order of A^,), a clear trend towards 
smaller intermittency is present. It is interesting to note 
the approximate collapse of the data from the two rota- 
tion rates, suggesting that is the relevant time scale 
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FIG. 4: Intermittency factors 7^ = g/2 — Cg/C2 (same data as in 
Fig. [sjl. Error bars are only shown for the last curve (*) for clarity. 
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FIG. 5: Time evolution of the 4th and 8th-order intermittency 
factors for the two experiments, o, Q = 1.13 rad s~^; □, Q = 
2.26 rad s~^. The dashed lines show the intermittency factors 
from the She-Leveque model, 74 = 0.161 and 73 = 0.824. The 
error bars (not shown) are of the order of the scatter, A74 = 0.05 
and A78 = 0.14. 



for the intermittency reduction. 



The fact that the factors 7^ start decreasing from the 
beginning of the decay is probably due to the low instan- 
taneous Rossby numbers when t ~ Iq (see Table |l]). A 
crossover between constant 7^ at early time and a de- 
crease at larger times would be actually expected for 
larger grid Rossby number Rog. However, a large grid 
Reynolds number Rcg is required for a developed turbu- 
lence to remain throughout the self-similar decay regime, 
up to the Ekman cutoff t tc, limiting the maximum ini- 
tial Rog at fixed rotation rate and grid size. We finally 
note that, extrapolating the trend towards 7^ — > in 
Fig. [5] suggests that the upper bound tc of the self-similar 
decay regime, in our experimental conditions, prevents 
from a clear observation of a vanishing intermittency, 
which may occur after 10-15 tank rotations. 

To summarize, our measurements of SF in decaying ro- 
tating turbulence show a strong increase of the exponents 
C,q during the decay, which essentially follows the increase 
of the second-order exponent C2 • It is worth noting that 
values for C2 larger than 1 are found, in contradiction with 
the S2{r) ^ r [i.e., E{k) ~ fc~^] phenomenological law 
for rotating turbulence, derived under the assumptio n of 
nonlinear interactions governed by the timescale f^-i 12115] 
Once normalized by C21 a marked increase of Cq/C2 is ob- 
served, a clear signature of a reduction of intermittency 
induced by the background rotation. This intermittency 
reduction is comparable to the one reported in a forced 
DNS by Miiller and Thiekjl^l but it is less pronounced 
than in the forced experiment by Baroud et aZ.j'^' proba- 
bly because of the limited temporal range of self-similar 
decay due to the Ekman dissipation regime, which is spe- 
cific to the decaying case. 
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